First-Principles Study of Integer Quantum Hall Transitions in Mesoscopic Samples 
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We perform first principles numerical simulations to investigate resistance fluctuations in meso- 
scopic samples, near the transition between consecutive Quantum Hall plateaus. We use six-terminal 
geometry and sample sizes similar to those of real devices. The Hall and longitudinal resistances 
extracted from the generalized Landauer formula reproduce all the experimental features uncovered 
recently. We then use a simple generalization of the Landauer-Biittiker model, based on the inter- 
play between tunneling and chiral currents - the co-existing mechanisms for transport - to explain 
the three distinct types of fluctuations observed, and identify the central region as the critical region. 
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Although the Integer Quantum Hall Effect (IQHE) 
is a generally well understood phenomenon, recent ex- 
periments on mesoscopic samples Q, Q uncovered unex- 
pected behavior in the seemingly noisy fluctuations of 
the Hall {Rh) and longitudinal [Rl) resistances. Pre- 
viously, fluctuations in resistance had been observed in 
mesoscopic samples with a phase coherence length 
larger than the sample size |^ Q |^ ; they are totally ran- 
dom, similar to universal conductance fluctuations 0. 
By contrast, Peled et. al. find P, Q that the transi- 
tion between the n*^ and (n -I- 1)"* plateaus of the IQHE 
has three distinct regimes: (i) on the high-_B side, both 
Rh and Rl have large but correlated fluctuations, such 
that i?i + Rii = h/ne^; (ii) for intermediate B values, 
Rh and Rl continue to exhibit fluctuations, but their 
sum is no longer constant; and (iii) on the low-_B side, 
Rh = h/{n + l)e^ has no fluctuations, whereas Rl still 
does. Moreover, Rl + Rh — R2t holds for all B val- 
ues 13]. Changing the sign of the magnetic field B ^ —B 
also has interesting consequences, as discussed later. For 
n = 0, regions (i) and (ii)are replaced by the transition 
to the insulating phase |l|. In this Letter, we explain 
the physics behind these observations in a unified theory, 
and analyze the implications for further experimental and 
theoretical study. 

The relation Rl + Rh = R2t was first proposed by 
Streda et. al. 0, while the fluctuations of regime (iii) 
are reminiscent of Jain and Kivelson's theory on the resis- 
tance fluctuations of narrow samples These theories 
were questioned by Biittiker based on formulas ap- 
propriate for a four-terminal geometry We take an 
approach similar to Biittiker's and use the multi-probe 
Landauer formula 0, ^| to calculate the resistances 
measured experimentally. However, we mirror the ex- 
periments by including six contacts in our model, namely 
the four voltage probes plus the source and the drain for 
the electrical current. Our model enables us to calcu- 
late both Rl and Rh, and reveals the very rich physics 
underlying the observation of the mesoscopic IQHE. 

The response function of the system is a 6 x 6 con- 



ductance matrix g, with which the current-voltage rela- 
tion reads la — J2i3 dap'^P- Here, la is the out-going 
current on the contact a = 1, • • • , 6 and Va is the corre- 
sponding voltage. For a mesoscopic Hall bar, g charac- 
terizes the electrical response (since the system is inho- 
mogeneous, local quantities such as the current density 
are not conceptually well-defined), g is calculated [llj 
by solving a multi-channel scattering problem: ga,f3^a = 
IT X^i j \'tai,fijf: where tai,i3j is the transmission ampli- 
tude from the j*^ transverse channel of contact /3 into the 
i*^ transverse channel of contact a for an electron at the 
Fermi energy Ep. Due to the absence of charge accumu- 
lation and to gauge invariance, J^adap = 12p9ap — 0- 
As a result, diagonal gaa are uniquely determined by 
gaa = - J2i3^a9a0 = " Y^f3^a9(3a- This imposcs a Con- 
straint on the off-diagonal elements of g for each a. 

Our model is sketched in Fig.^a). Six perfectly con- 
ducting contacts are linked to a 4/imx2/im sample which 
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FIG. 1: (color online) (a) Typical potential Vd -I- 14 of a 
4/imx2/im sample, (b) Chiral (arrows) and tunneling (resis- 
tors) currents in our model. This direction of chiral currents 
corresponds to B entering the page, (c) Some semi-classical 
current distributions parameterized by our model, (d) Jain- 
Kivelson tunneling for high-i^. See text for further details. 



has a disorder potential Vd(r) and a background potential 
Vf,(r). Vd{r) is a sum of random Gaussian scatterers gen- 
erating elastic scattering inside the sample, while Vb{r) 
confines the electrons to the sample and creates edge 
states. In our simulations, we include LxLyB/(j)o ~ 10^ 
states of the lowest Landau level (LLL), where L^Ly 
is the area of the sample and 0o = h/e is the mag- 
netic flux quantum. The sample Hamiltonian is a large, 
sparse matrix Hnm — {ipnlVb + Vd\ipm) ■ Contacts are 
modeled by ensembles of one-dimensional tight-binding 
chains attached to localized eigenstates on the corre- 
sponding edges of the sample. We have verified the multi- 
probe Landauer formula for our model using the linear 
response theory. This derivation and further modeling 
details will be reported elsewhere For a given mag- 
netic field B, we numerically solve the scattering problem 
for different values of the Fermi energy and obtain g. The 
filling factor v is also a function of Ep, and thus we can 
find the dependence of the conductance matrix g on v. 

The resistances are then computed from g. In the usual 
setup the current is injected in the source and extracted 
in the drain -h = h = I; hi = (-/ / O)^. 
Without loss of generality we set 1=1 and 14 = 0. 
The other five contact voltages are uniquely determined 
from /i4 ~ g -V . We define two longitudinal resistances 

^f4,23 = (^2 - V^)/I = V2- ^3, i?f4,65 = ^6 " ^5, and 

two Hall resistances i?i4_62 = ^ — ^2, ^14,53 = V5 — V3. 

In Fig. 12 we plot representative matrix elements gap 
as a function of v. For v > 0.5, we find ga,a+i — * e'^/h 
(ifa = 6,a-|-l = l), with all other off-diagonal matrix 
elements vanishing. In other words, all electrons leaving 
contact a + 1 arrive at contact a. It follows that here 

ga/3 5q/3 = -J^ + + ^aG^fil) ■ (1) 

iFrom /i4 = g^°^-V we find V5 = Ve = h/e^, V2=V3 = 0, 
thus i?f4,62 = i?f4,53 = h/e^ R^,^,, = R^,^,, = 0. This 
shows that the first quantized plateau is due to the chiral 
edge currents [shown as oriented thick lines in Fig.^b)], 
which become established for v > 0.5. Variations of g{i') 
from g^°^ give rise to fluctuations in the resistances. From 
Fig. 13 we also see that if i/ < i/^ (indicated by the vertical 
line), gafi ~ 5/3q with high accuracy, i. e. g is symmetric. 
For v > i^c, g is no longer symmetric. The reasons for 
this behavior and its consequences are discussed later. 

Using the conductance matrix g{i') plotted in Fig. 13 
we now compute the values of the various resistances as 
a function of < < 1. Fig. shows a pair Rl and 
Rh- Three different regimes are clearly seen: for v > 
0.46, Rh = h/e^ and Rl = 0, corresponding to the first 
IQHE plateau. For 0.42 < v < 0.46, Rl exhibits large 
fluctuations, however Rh is still well quantized. This is 
precisely the type of behavior observed in Ref. Q. For 
V < 0.42 the transition to the insulating phase occurs, 
and both resistances increase sharply. The fluctuations 




FIG. 2: (color online) Representative conductance matrix el- 
ements, in units of /h, as a function of the filling factor u. 
The left (right) panel shows 523,545 and 552, respectively 532, 
554 and 526 characterizing transport in the (against the) di- 
rection of the edge currents. Results are almost identical on 
the left of the dot-dashed line, but different on its right. 



are very large and narrow because the calculation is done 
at T = 0. At finite T, the peaks are smeared out. 

The transition 1 < ly < 2 can also be simulated us- 
ing the same (/{v) matrix of the LLL. Similar to Ref. 
we assume that the completely filled spin-up LLL con- 
tributes its background chiral edge current. As a result, 
we simply add 5*^°^ = 5(1/ = 1) of Eq. to the values 
of g(i^) describing the partially filled spin-down LLL. Al- 
though the two LLLs have different spins, the contacts 
mix electrons with both spins in equilibrium, justifying 
this addition. Resistivities Ru:.e2 and i?f4 23 computed 
for 5^°) -I- g{L') are shown in panel (b) of Fig. |31 whereas 
in panel (c) we plot their sum. The three regimes found 
experimentally 0, are clearly observed. At low-v 
(high-B), the fluctuations of Rh and Rl are correlated, 
Rl + Rh = h/e'^. At high-j/ (low-B) Rh = h/2e^ is 
quantized while Rl still exhibits strong fiuctuations. In 
the intermediate regime, both Rh and Rl have strong, 
uncorrelated fluctuations. The other pair, R^^ 53 and 
i?f4 g5, also exhibits the three regimes, although their de- 
tailed fluctuations are different from -Rf^ g2 and i?f4 23- 
From over 20 different simulations we found that the low- 
V regime where Rl + Rh = h/e^ is a very robust feature, 
although it is maintained up to different values of v in 
different samples. The high-j/ regime with fluctuations 
in Rl and quantized Rh is seen frequently. However, 
when strong direct tunneling occurs between the source 
or the drain and their nearby voltage probes, Rh also 
fluctuates. Such strong tunneling is an artifact of our 
simulation 0. We suppress it by isolating the source 
and drain from nearby contacts with triangular poten- 
tial barriers in the corners of the sample [see Fig. ^a)]. 
Figure ISIc) also compares i?{4^23 + ^14,62 = -Ri4,63 with 
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FIG. 3: (color online) Rl and Rh calculated from the con- 
ductance matrix displayed in Fig. |21 in units of h/e^. (a) 
Transition from the insulator to the first IQHE plateau in 
the LLL. (b) Transition from the first to the second IQHE 
plateaus, (c) The sum Rl + Rh of the resistances shown in 
(b), and R2t — 0.5/i/e^. Vertical lines indicate the boundaries 
of the critical region. See text for further details. 

i?2t = i?63.63- [In the setup for measuring i?2t, the cur- 
rent is /63 = (0 1 -1)^, and i?2t = V^e - V3]. As 
found experimentally 0, the two curves are very similar. 

So far, we have demonstrated that our numerical sim- 
ulations recapture faithfully the experimental results. 
We now explain the underlying physics using a sim- 
ple but very general model. For the given constraints 
and using logical induction, we find 12] that g can be 
decomposed as a sum over loops connecting contacts: 
g = ^c(fci,--- , A:„)f(fci, • ■ • ,kn). Here, c are positive 
numbers and • • • , fc„) — l{ki, + • • • + l{kn, fci), 
where la(3{a,b) = ^ {SaaS/3b ~ ^S^aS/sa - ^SabS/3b) con- 
tributes to a single off-diagonal element gab- A two- vertex 
loop r(a,6) describes ah/e^ resistor between contacts a 
and b. Since rap{a,b) — rfja{a,b), these terms are the 
symmetric part to g. The asymmetric part of g describes 
chiral currents, whose direction of flow is dictated by the 
sign of B. In particular, f(l, 2, 3, 4, 5, 6) — of Eq. 
n describes the edge currents of a LL, but shorter chiral 
circuits may also develop at intermediate fillings v. 

At low-j/, all states are localized and transport in the 
LL can only occur through tunneling. Consider the semi- 
classical case sketched on the left side of Fig.QIc). Elec- 
trons can go from 2 to 1 either through direct tunneling 
(probability ti2), or they can tunnel to the localized state 
near contact 6 and from there back into 1, with proba- 
bility (1 — ii2)i26^i6(l ~ tu)- Electrons can make any 
number of loops before entering 1, the total sum being 

P2^l = ^912 = [tl2 + ^26*16 — 2ti2t26il6]/(l — ^12^16^26)- 

Similar arguments give pi^2 — ^921 — \t\2 - ^26^16 - 
^I2ii6]/(1 -^12*16*26)- 516,561,526 and 362 are computed 
similarly. We define Tab = inin(Pa^b,Pb^a) > 0, and 



Cab = max(pQ^;,,pb^Q) - Tab > 0. We find C12 = C26 = 
cgi = c = ti2ti(i + ti2t26 + ti6t26 + 0{t^); and up to O(t^), 
7*12 ~ ^12, '^16 ~ ^16 and r26 ~ ^26- These processes 
contribute a total of r-i2f(l, 2) + rigf(l, 6) + r2Gf{2, 6) -I- 
cf(l, 2, 6) to g. The symmetric resistance terms, of order 
t, are due to direct tunneling between contacts, and at 
low-iy they dominate the small chiral current, of order t^. 
This explains why for v < Vc, g symmetric with small 
off-diagonal components (see Fig. EJ. At high enough i^, 
edge states connecting consecutive contacts appear. As 
already discussed, as J/ — > 1, 5 ^ 5'°' = f(l, 2, 3, 4, 5, 6). 
For intermediate v, shorter chiral loops containing edge 
states can be established through tunneling, as sketched 
on the right side of Fig. ^c). Assume that an electron 
leaving contact 3 can tunnel with probabilities and 
to and out of a localized state, to join the opposite edge 
current and enter 5. It follows that ^3^5 = ^553 = 
^3^5/[l ^ (1 ^ ^3)(1 ^ ^5)], whereas P5_>3 = (no electron 
leaving 5 enters 3). Then rss = and the contribution 
to g is just ^3^5/(5, 3). This term combines with parts 
of i(3,4) and /(4,5) to create a chiral current cf(3,4, 5), 
where c ~ Ps^s- Physically, this represents the backscat- 
tered current of the Jain-Kivelson model |^ . 

In general, one has to sum over many types of compet- 
ing processes, involving both tunneling and edge states, 
but g can always be decomposed into symmetric re- 
sistances plus chiral loops. Consider the general form 
g = ng(0)+ri2f(l,2) + ri6f(l,6)+r26r(2,6)+r34r(3,4) + 
r45r(4, 5)+r35r(3, 5)+co5(°)+cir(l, 2, 6)+C2r(2, 3, 5, 6) + 
C3f(3, 4, 5) + C4f(l, 2, 3, 5, 6) + C5r(2, 3, 4, 5, 6). The first 
term describes the contribution of the n completely filled 
lower LLs. All other terms describe transport in the LL 
hosting Ef [see Fig. ^b)] , with the restriction that there 
is no tunneling between the left and right sides of the 
sample. This is justified physically because tunneling be- 
tween contacts far apart is vanishingly small. The largest 
such terms, 7-23 and r^Q, are found to be less than lO"'^ 
[see e.g. Fig. El, where r23 = h/e^ ■ min(532, 532)]. We 
solve both /14 = g-V and Iqs = g-V and find the identity 

^14,63 — ^63,63 — ^ — ; ; ; 

n + Co + C2 + C4 + C5 

Since i?63,63 = R2t, whereas i?i4,63 = ^14,62 + ^14,23 = 
Ri4,Q5 + -R14.53, this means that i?2t = Rh + Rl irre- 
spective of the value of the 12 parameters. In other 
words, this identity is obeyed for all v, in agreement with 
Fig. inic) . (Adding r23 and terms leads to perturba- 
tionally small corrections jl2i|'). Here n+co + 02 + C4 + C5 is 
the total chiral current flowing along the 6 — > 5 and 3 — > 2 
edges. In particular, at low-J^ chiral currents in the LL 
hosting Ep are vanishingly small cq — C2 = c^ = c^ = 
(there are no edge states established yet, and pure tun- 
neling contributions are of order t"^, as already discussed. 
Below Vc, all t < 10"'*, see Fig. It follows that here 
Rl + Rh = h/ne^ , explaining the perfect correlations in 
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the pattern fluctuations at low-r^ of the two resistances, 
observed both experimentally and numerically. 

The high-!/ regime with quantized Rh and fluctuating 
Rl can also be understood easily. As already discussed, 
transport in the LL hosting Ep is dominated here by 
the edge states; tunneling between opposite edge states 
(facilitated by localized states inside the sample) creates 
backscattered currents, as in the Jain-Kivelson model 
We sketch this situation in Fig. ^d), where ti, t2 re- 
spectively include all possible tunneling processes lead- 
ing to backscattering on the corresponding pairs of edge 
states. Reading the various probabilities off Fig.^d), we 
find that g = ng^") + (1 - <i - 12 - is)^^"^ + h[rih 2, 6) + 
f(3,4,5)] + t3''(l,2,3,5,6) + iif(2, 3, 4, 5, 6). The first 
term represents the contribution of the lower n com- 
pletely filled LLs, the others are the forward and the 
backscattered chiral currents in the LL hosting Ep. /14 — 
g-V is trivial to solve. We find R14Q2 — ^14.53 = (^* + 
l)e^, i. e. the Hall resistances are precisely quantized, 
whereas R14 23 = ^14,65 = [^/("- + • h/in + 1- ^2). 
Since t2 has a strong resonant dependence on Ep (or v), 
it follows that Rl fluctuates strongly. In particular, if 
n = (transition inside spin-up LLL), R^ can be ar- 
bitrarily large when t2 ^ 1, whereas in higher LLs the 
amplitude of fluctuations in Rl is h/{n+ l)e^ or less, as 
observed both experimentally and in our simulations. 

If B changes sign, we have verified that the identity 
g{-B) = holds [uj. The reason is that time- 

reversal symmetric tunneling is not affected by this sign 
change, while the flow of the chiral currents is reversed. 
The model mirrors itself with respect to the horizon- 
tal axis if B changes sign, see Fig. ^ The solutions 
of 7i4 = g{—B) ■ V are then related to the solutions of 
/i4 = g{B) ■Vhyv2 = Vq, V3 = V5, V5 = V3 and vq = V2, 
provided that the same index exchanges 2 <-> 6, 3 ^ 5, 
are performed for all Vab- Terms not invariant under 
this transformation are ri2, rig, r43, r45, r23 and rsg. 
As already discussed, the last two terms are vanishingly 
small. In the experimental setup, the first four terms are 
also very small, due to the long distance between source 
and drain, and their nearby contacts jl5l |. The dom- 
inant terms r26 and ras are invariant under the index 
exchange. It follows then that Ri4^23{^) = ^u.esl"^) 
and vice versa, i. e. with good accuracy, the fluctuation 
pattern of one Rl mirrors that of the other Rl when 
B — > —B. This symmetry has indeed been observed ex- 
perimentally, with small violations at low-i/ due to 
perturbative corrections from the non-invariant tunnel- 
ing contributions ri2 — rig and r43 — r45. 

We now summarize our understanding of the various 
results of IQHE measurements on mesoscopic samples. 
Similar to experiments, we find that the transition in 
higher LLs is naturally divided in three regimes. At low- 
v, the LL hosting Ep is insulating and there are no edge 
states connecting the left and right sides of the sample. 
If tunneling between left and right sides is also small. 



we find that the fluctuations of pairs of resistances are 
correlated with excellent accuracy, Rh + Rl = h/ne^. 
This condition is obeyed if the typical size of the wave- 
function (localization length) is less than the distance 
between contacts 2 and 3. When the localization length 
becomes comparable to this distance, an edge state is es- 
tablished and the correlation between Rl and Rh is lost. 
On the high-!/ side the edge states are established, but 
localized states inside the sample can help electrons tun- 
nel between opposite edges, leading to back-scattering 
like the Jain-Kivelson model. In this case, we showed 
that Rl fluctuates while Rh is quantized. Tunneling be- 
tween opposite edges is likely only if the typical size of 
the wave-functions is slightly shorter than the distance 
between opposite edges. It is then apparent that the 
central regime in Figs. 121(b) and (c) corresponds to the 
so-called "critical region", where the typical size of the 
electron wave- function is larger than the sample size (dis- 
tance between contacts 2 and 3, at low-i/, or between 2 
and 6 at high-i/). In these mesoscopic samples, the volt- 
age probes act as markers on a ruler, measuring the local- 
ization length of the wave-functions at the Fermi energy. 
To our knowledge, this is the first time when the bound- 
aries of the critical region are pinpointed experimentally. 
This opens up exciting possibilities for experimentally 
testing the predictions of the localization theory. 

To conclude, we used both first-principles simulations 
and a simple model to explain the phenomenology of the 
mesoscopic IQHE, for six-terminal geometry. We iden- 
tified tunneling and chiral currents as coexisting mech- 
anisms for charge transport in mesoscopic samples, and 
argued that the boundaries between the three distinct 
regimes mark the critical region. 
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